Chiral molecular films as electron polarizers and polarization modulators 
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Recent experiments on electron scattering through molecular films have shown that chiral 
molecules can be efflcient sources of polarized electrons even in the absence of heavy nuclei as 
source of a strong spin-orbit interaction. We show that self-assembled monolayers (SAMs) of chiral 
molecules are strong electron polarizers due to the high density effect of the monolayers and explic- 
itly compute the scattering amplitude off a helical molecular model of carbon atoms. Longitudinal 
polarization is shown to be the signature of chiral scattering. For elastic scattering, we find that at 
least double scattering events must take place for longitudinal polarization to arise. We predict en- 
ergy windows for strong polarization, determined by the energy dependences of spin-orbit strength 
and multiple scattering probability. An incoherent mechanism for polarization amplification is pro- 
posed, that increases the polarization linearly with the number of helix turns, consistent with recent 
experiments on DNA SAMs. 



The surprising experimental finding that self assem- 
bled monolayers (SAMs) of chiral organic molecules on a 
gold substrate can act as highly efficient spin filters or po- 
larizers, with a polarization factor larger than 60% dl, has 
prompted us to examine the underlying physical mech- 
anism responsible for this behavior. In the absence of 
external fields, spin-dependent effects in electronic trans- 
port across an interface between two materials are es- 
sentially determined by electron-electron correlation and 
spin orbit (SO) interaction!^. In the experiments we are 
referring to, photoelectrons emitted from the gold sur- 
face are transmitted through the SAM at energies larger 
than the interfacial tunneling barriers [H Under these 
conditions electron transmission is mostly a one-particle 
process and therefore the explanation for the spin polar- 
ization effects must physically depend on either the pres- 
ence of strong internal magnetic fields or the spin-orbit 
interaction. We argue in this letter that the most strik- 
ing experimental observations: the role of chirality and 
the prevalence of longitudinal polarization, can be under- 
stood from a basic principles model using a one particle 
hamiltonian that includes the bare spin-orbit interaction. 

Spin-dependent electron scattering has been observed 
in gas phase experiments Q with chiral molecules as tar- 
gets, and the observed asymmetry factor for spin polar- 
ization can be well understood as due to the result of 
a substantial spin-orbit interaction caused by a heavy 
atom present in the molecule. A similar explanation can 
be attempted for the experiments with SAMsQ, but the 
objection has been raised that for the low atomic num- 
ber atoms, C, H, N and O present in the monolayers, 
the magnitude of the SO interaction is too small to ex- 
plain the observed polarization[6|. In a very recent exper- 
iment reported by Gohler et al[l|, a direct measurement 




FIG. 1: Scattering model of a chiral assembly of molecules 
chemisorbed on a gold surface. Chiral molecules are assumed 
helices of radius a and pitch b built from carbon atoms, whose 
electronic state is fully described in the model. The electrons 
are incident in the z direction and are assumed to be scat- 
tered internally within the layer, and then received at a Mott 
detector. Electron wavefunctions from surface are normalized 
within a rectangular box of volume v. 

of spin polarization with a strong longitudinal polariza- 
tion was obtained for an unpolarized emerging electron 
beam. This experiment establishes in a definitive way 
that the SAM, composed of light atoms only, can polar- 
ize electrons regardless of the initial origin and magni- 
tude of the electron polarization. We use a model for 
elastic spin-dependent multiple scattering through chiral 
molecules to explain these experimental findings. Three 
key ingredients emerge: (i) Although transverse polar- 
ization can arise under inelastic events from point tar- 
gets, the appearance of longitudinal polarization requires 
the existence of chirality. (ii) Strongly enhanced spin 
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orbit interaction is due to the fact that the density of 
molecules in a SAM is comparable to that of a solid 
structure, thereby increasing the electron wavefunction 
overlap with deeply penetrating orbitals as compared to 
the gas phase (hi) Multiple scattering is needed to 
produce longitudinal polarization, so that there is an en- 
ergy range where electron polarization is most effective in 
the few electron volts range, as seen in the experiments. 
We claim that sources of spin-orbit coupling analysed^ 
such as electric polarization fields, even though as large as 
10^ — 10''V/cm[8|), cannot account for the spin-splitting 
energy necessary for strong spin polarization. Such of 
SO also do not carry the field configuration to transfer 
the appropriate longitudinal angular momentum to the 
transmitted electrons. It is the wavefunction penetration 
to the nucleus in atomic orbitals that will yield electric 
fields of 10^°V/cm at distances of a Bohr orbit within the 
atom. This is the source of the peculiar strength of SO 
in non-centrosymmetric semiconductors 

The detailed calculation of the scattering cross sec- 
tion and polarization of even a small molecule like H2O 
[loj can be very demanding, taking into account static, 
exchange and spin-orbit interactions in detail. In or- 
der to capture the coherence effects by connecting them 
to molecular parameters, it is necessary to have ana- 
lytical solutions to simpler models. In this paper we 
solve explicitly for the scattering of electrons from a sin- 
gle oriented helical moleculeQ, in a SAM as depicted 
in Figll] in the Born double scattering approximation. 
We consider the elastic scattering of spin-quantized elec- 
trons from a chiral potential with explicit consideration 
for the SO interaction and scalar potential effects from 
the particular atomic centers that conform the molecu- 
lar structure. Wc solve for the transition amplitude in 
the scattering formulation where electrons arc incident 
with relative momentum fc, and arc subject to the po- 
tential V{r, cr) ~ V{r)l + a&. [W{r) A p], at position r, 
including both scalar and spin dependent contributions 
with a ~ — (2ni^e)^ hare Pauli SO prefactor and V{r) 
the atomic Coulomb potential. 

In order to compute the scattering amplitudes, we need 
to expand the t matrix in powers of the potentials in- 
volved. In the double scattering approximation, and as- 
suming non-overlapping multicentral potentials 11[ , we 



have that for a multiatomic molecule assumed to be com- 
posed of carbon atoms C*-^-* through C*^^^ 
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which includes both single and double scattering con- 
tributions. Here, k^, k;, are the scattering incom- 
ing/outgoing wavevectors. The subindices of the single 



scattering t matrix correspond to the particular atom in- 
volved C^'^ and Go = exp{ik\rA — rB|)/47r|r^ — rsj for 
atomic centers A and B. The single scattering tA matrix 
elements are given by (ko|t^|ki) = exp(i(ki - ko) • 

ryi)fyi(ki, ko) with the single scattering amplitude 
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where ko, k^ stand for collision input and output wa vec- 
tors between any pair of atoms. The potentials have been 
translated to the origin, as their phase factors have been 
extracted from the integral. The denominator v repre- 
sent a volume per molecule where the incident electron 
wavefunction is normalized, and T{ki,ko) is the expec- 
tation value of the atomic potential in the plane-wave 
basis. 

Note a strong result from the above expression; no lon- 
gitudinal polarization can result from single scattering 
events no matter the symmetry of the potential since only 
transverse components can arise (Pauli matrix compo- 
nents in the scattering plane). Nevertheless, double scat- 
tering events that result from products of single scatter- 
ing amplitudes will contain all polarizations, as a result 
of the algebra of the Pauli matrices, physically control- 
ling the transfer of orbital to spin angular momentum. 
Thus, one needs to go beyond single scattering to gain 
longitudinal polarization. Higher order scattering events, 
as we will see, have energy requirements that should be 
reflected in polarization experiments and will serve to 
prove their role. We also argue that higher order scatter- 
ing is the only way to get chirality imprinted on the lon- 
gitudinal spin signature. In view of the previous points, 
we consider the molecular potentials as a superposition 
of their atomic nuclei, disregarding bonding effects 
Wc employ the wavcfunctions of Clementi and Roetti 



[12| 

131 



derived from a Hartrcc-Fock (HF) model, as convenient 
expressions to simplify explicit calculations. 

The situation of highly packed SAMs can be compared 
to that of semiconductors, in the sense that the major 
part of the SO interaction strength comes from prob- 
ing electric fields close to the atomic cores. The SO re- 
tains the strength associated with individual atoms in 
the solid 0, the spin splitting energy being very similar 
(0.006 eV for carbon[l^). We can correctly estimate the 
SO interaction, adjusting no parameters i.e. from the mi- 
croscopic Hamiltonian, by calculating the single particle 
amplitude as follows: We normalize the electron incident 
plane wave in a box whose volume is equivalent to the vol- 
ume per molecule in the SAM. We then expand the plane 
waves in the local atomic wavefunction basis; T{ki,ko) = 

J2nl,n'l'i^o\Vnl){'Pnl\V{r)\ipn'l'){V>n'l'\ki), with \ipnl) the 

HF wavcfunctions and |fci,o) the input-output normal- 
ized plane wave vectors, and V{r) the bare Coulomb po- 
tential. This estimate has to be done for each incom- 
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ing/outgoing wavevector combination, to properly ac- 
count for momentum transfer and adding over the avail- 
able orbitals of the target atoms which for carbon in- 
clude 2p and higher levels with I ^ that can contribute 
to the SO scattering strength. This procedure is a very 
simple and robust approach that accounts for the correct 
spin splitting and the density effect of the scattering tar- 
gets, and for the order of magnitude differences between 
gas phase polarization^] and that resulting from electron 
scattering in the SAMs. 

It is useful to rewrite the full scattering amplitude (one 
and two body scattering Eq. [Ij in the form 
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where the functions ft.j;(ka, k;,) are to be determined. The 
scattered polarization component Pj, when the incident 
state is a completely unpolarized state p = (1/2)1, can 
be computed from the output density matrix p' = f pf ^ 
through the relation Tr(crp')/Trp'. In terms of the ex- 
pansion in Eq|31 the polarization can then be written as 
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where the indices (/, m, n) are taken in cyclic order and Jft 
denotes the real part of the argument. We immediately 
note a central feature of polarization effects in electron 
scattering: It is an interference effect between the spin- 
orbit interaction hi with i ^ and the purely potential 
term ho. One needs the presence of both to guarantee 
polarization. 

In Figl2]we depict the electron polarization as a func- 
tion of the pitch b of the discrete helix built from six 
carbon atoms, that comprise a single turn of the struc- 
ture. The range of k values considered is between and 
1 a.u. (energies between and 10 eV), in the range that 
has been explored in polarization experiments [l|. For 
all values of the incident fc, the longitudinal polarization 
vanishes as chirality, as measure by the pitch b, goes to 
zero. The relation between longitudinal polarization and 
pitch constitutes a measure of chirality itself, although 
the measure is obviously non-monotone, since there are 
important interference effects at play. The longitudinal 
polarization, can be of any sign depending on the range 
of wavevectors considered [1] and the magnitude of the 
pitch of the scattering target. Note that one can achieve 
polarizations of 1% for k = 1 a.u. which is three orders of 
magnitude larger than gas phase polarizations reported 
with heavy atoms and compatible with a comparable 
molecule in a SAM, such as Lysine A molecule such 
as DNA contains a much larger number of atoms per 
turn and a double helix. It is beyond the scope of our 
approach to include such details in our model. 
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FIG. 2: Total longitudinal polarization per molecule, aver- 
aged over outgoing scattering angles, versus helix pitch b as 
defined in the inset, for representative values of k. No longitu- 
dinal polarization is obtained for b = (chirality zero) for any 
value of the wavevector. The polarization is a non-monotonic 
function of b beyond a certain pitch. The inset (bottom panel) 
shows the change of the sign of the longitudinal polarization 
when chirality changes sign. 



To assess the energy dependence of the polarization, 
the best measure is the scattering probability in terms of 
the components defined by Eql31 since polarizations and 
asymmetries always yield a relative quantity to the total 
scattered intensity. In Figl3]we see |/i3p i.e. the ampli- 
tude of the longitudinal polarization, as a function of the 
incident wavevector k and scattering angle 0. One can 
readily sec that the relevant energy regime for chirality 
effects, in the double scattering approximation, are in the 
units of eV range (in a.u. in the figure). Since longitu- 
dinal polarization requires at least double scattering, the 
electrons must be deflected at large angles for intermedi- 
ate dispersion events. Nevertheless, large deflections can 
only occur with the denominator handicap oc jfco — fci|^ 
the momentum transfer. If k is increased by a factor of 
two, the spin-orbit interaction effects can be reduced by 
an order of magnitude, concomitantly reducing longitudi- 
nal polarization. There is a corresponding effect at lower 
k where the spin-orbit interaction decreases linearly. So 
there will be a peak range in energies where longitudinal 
polarization is observed. 

Up to this point we have only computed the polarizing 
power of a single turn of the helix under double scatter- 
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FIG. 3: The unnormalized probability amplitude for longitu- 
dinal polarization as a function of the incident wavevector and 
angular dependence. The figure shows the window of wavec- 
tors (in figure between k=0.3 and 1 a.u.) where longitudinal 
polarization is observed. 



ing conditions. Successive turns of tlie Irclix or additional 
layers of cliiral molecules can be treated incoherently by 
assuming layers of one helix turn thick. Since all scat- 
tering occurs at room temperature in the experiments, 
one does not expect long coherence lengths beyond a few 
nanometers. The scattering is then presumed to proceed 
by successive incoherent double scatterings through dif- 
ferent turns of the helix. This can involve either the same 
helix or the ones within the lateral distance defined by 
the transverse component of the wavevector of the elec- 
trons emitted from the surface. 

If the polarizing power of a single turn of the helix is 
Punpoi = Tr[crfft]/Tr[fft] i.e. that given by Eqgl then 
the polarization through an additional turn of the helix 
(if the turns add up incoherently) , can be computed from 
the density matrix = f(l + P(*) • (T)ft/2. Ignoring 

angular dependencies of incident electrons from one helix 
turn to the next, one obtains 
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where P^*-* is the polarization achieved after the z-th turn 
of the helix and h = (/ii, /12, /i3)(see Eq[3l). Figure |4] 
shows the polarization as a function of the number of 
helix turns for two different pitches smaller than the he- 
lix radius. We see a linear growth of the longitudinal 
polarization, proportional to the pitch size. The linear 
behavior, in spite of the complex relation between differ- 
ent components of the polarization, lends support for in- 
coherent amplification. For larger pitch values, however, 
the relation between different components of the polariza- 
tion can make the amplification non-monotonous, chang- 
ing the polarization sign as the number of turns increase. 
The same mechanism can be envisaged for stacked lysine 
SAMsil. 

Summarizing, we have presented an elastic scattering 
scenario to describe spin-orbit interaction based polariza- 
tion of electrons in chiral SAMs. The crux of our model 
is to consider multiple scattering off a helical model built 
from carbon atoms, disregarding bonding effects. The 
analytical approach, albeit consideration of only dou- 
ble scattering events, allows understanding the under- 
pinnings of longitudinal polarization, the energy windows 
involved and the role of the chiral parameter. We have 
accounted for the strength of the SO interaction from 
the microscopic interactions by projecting the normal- 
ized plane waves onto the available atomic orbitals of 
the carbon scattering centers. This approach has no ad- 
justable parameters and accounts directly for the SAM 
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FIG. 4: The longitudinal polarization as a function of the 
number of turns in the helix for b < a — 1.88 a.u. and k = 
1 a.u. The polarization grows linearly with helix length, and 
a pitch dependent slope. The inset depicts the incoherent 
scattering from separate slabs of one turn in thickness. 



density. Double and higher order scattering will be sig- 
nificant only within an energy window due to momentum 
transfer denominators that hinder the lateral scattering 
necessary to explore the target molecule. Finally, the 
short coherence lengths for electrons at room tempera- 
ture, lead to the scenario that the linear increase of the 
longitudinal polarization with chiral chain lengths [l| is 
due to incoherent multiple forward scattering. 
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